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Abstract 

We analyze the structure of the Coulomb and unitarity corrections to the single e~^e~ pair 
production as well as the cross section cr„ for the multiple pair production in collision of ultra- 
relativistic nuclei. In the external field approximation we consider the probability of one pair 
production at fixed impact parameter p between colliding nuclei. We obtain the analytical result 
for this probability at large p as compared to the electron Compton wavelength. The energy 
dependence of this probability as well as that of o"„ differ essentially from widely cited but incorrect 
results. We estimate also the unitary corrections to the total cross section of the process. 
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Recently the process of e~^e~ pair production in ultra-relativistic heavy-ion colhsions 
was discussed in numerous papers. This is connected with the beginning of operation of 
Relativistic Heavy Ion CoUider (RHIC) with the Loretz factor 7 = 108 and charge number 
of nuclei Z = 79. New collider LHC is scheduled to be in operation in the nearest future, with 
7 = 3000 and Z = 82. The cross section of one pair production in the Born approximation 
was obtained many years ago |^ and reads 

cTBorn = ^ flv^ - 2.198 + 3.821 L - 1.6321 , (1) 

27tt ni'^ ^ J 

where 

C = {ZAa)\ZBaf , L = H^aIb) ■ (2) 

Here a is the fine-structure constant , m is the electron mass, Z^ b are the charge numbers of 
the nuclei A and B. The nuclei A and B are assumed to move in the positive and negative 
directions of the z axis, respectively, and have the Lorentz factors '~fA,B- For the sake of 
simplicity we consider the process in the frame where both nuclei have the same Lorentz 
factor 7a = 7b = 7- In the present paper we consider the case L ^ 1. Since cross section (|l]) 
is huge, a pair production can be a serious background for many experiments. Besides, it is 
also important for the problem of beam lifetime and luminosity of colliders (see review 0). 
It means that the various corrections to the Born cross section, as well as the cross section 
cr„ for n-pair production, are of great importance. At present, there is a lot of controversy 
in papers devoted to this subject (the corresponding references and critical remarks can be 
found in ^, P). Although some of the corrections were obtained correctly in some special 
cases, the whole consistent picture of pair production is absent. In the present paper we 
intend to elucidate some points in this problem. 

For 7 ^ 1 it is possible to treat the nuclei as sources of the external field, and calculate the 
probability of n-pair production Pn{p) in collision of two nuclei at a fixed impact parameter p. 
The corresponding cross section cr„ is obtained by the integration over the impact parameter: 

an= [ (fpPn{p). (3) 



The quantity which is important in further consideration is an average number of the pro- 
duced pairs at a given p: 

oo 

W{p) = Y.nPr,{p). (4) 
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The closed expression for W{p) was obtained in 
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Here p and ep (g and e^) are the momentum and energy of the electron (positron), u{p) 
and u{—q) are positive- and negative-energy Dirac spinors, ol = 7^7, 7± = 7° ± 7^, 7^ are 
the Dirac matrices, p± = tp ±p^, g-|- = ± g^, fc is a two-dimensional vector lying in the 
xy plane, k' = p± — k, and the function J-'{k) is proportional to the electron eikonal 
scattering amplitude in the Coulomb field. 
The function W{p) defines the cross section 
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which is called "the inclusive cross section" in |jTD|. Let us stress that the usual definition 
of the inclusive cross section is different: 
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To obtain cr-r it is necessary to perform the regularization of the expression for W{p). One 
of the possible correct regularizations is given in [^] and reads 

^A,B{k) = 271" y dppJoipk) {ex^p[2iZA,BaKo{pa±)] - 1} , a± = (p± + g±)/(27) , (8) 

where Jo is the Bessel function and Kq is the modified Bessel function of the third kind. 
The Born cross section can be obtained by the replacement 



(9) 



k^ + al ' 

where T^{k) is the first term in the expansion of J^{k) in Za. After the regularization the 
cross section ctt can be presented in the form: 



c cc 

ax = Cr Born + 0"^ + (Tj^ , 



(10) 



where a^^ is the Coulomb corrections with respect to one of the nuclei, and a^'" is the 
Coulomb corrections with respect to both nuclei. In the main logarithmic approximation 



3 



they were obtained in 



= -^^L'[fiZAa) + fiZBa)] 



^t"" = -\Lf{ZAa)f{ZBa), (11) 
97r 

where 



oo -| 

/(x) = V — ^ . (12) 

The expression for ay coincides with that obtained in by means of Weizsacker-Wilhams 
approximation. The accuracy of the expression (|iy), (|llD is determined by the relative order 
of the omitted terms ~ {Z^ b'^Y/L'^- This accuracy is better than 0.4% for the RHIC and 
LHC colhders. 



In a set of pubhcations [O, 12, 13, 14| it was argued that the factorization of the multiple 



pair production probability is valid with a good accuracy 

Pn{p) = ^e-^(-) . (13) 



The factor exp{—W) is nothing but the vacuum-to-vacuum transition probability Pq = 
1— Pn- Strictly speaking, such a factorization does not take place due to the interference 
between the diagrams corresponding to the permutation of electron (or positron) lines (see, 
e.g., JlOl)- Nevertheless, one can show that this interference gives the contribution which 
contains at least one power of L less than that of the amplitude squared. Therefore, in the 
leading logarithmic approximation we can use the expression ([T3|). 

Let us represent the cross section ai of one pair production as follows 

= (XT + a^nit = J d^pW{p) - j d^pW{p) (l - e-^('')) (14) 

Thus, the difference between cji and ctt is due to the unitarity correction o"unit- The main 
contribution to the first term (cr^) comes from p ^ 1/m. As for the second term, a^^iti 
the main contribution to it comes from p 1/m. If Z^^bo. ~ 1, then the exact W{p) at 
p ~ 1/m differs essentially from the Born result Wq{p). Since ar is studied in detail, here 
we investigate the unitarity correction and an for n>2. Below we consider two interesting 
cases: 

(i) Za,boi < 1, < 1 ; 

(ii) ZA,Ba<^lXL^l. 



Let us consider in detail the case (i), where it is possible to use Wq (obtained in the 
Born approximation) instead of the exact function W, expand the exponent exp{—W) in 
the unitarity correction and omit it in cr„ for n > 2. For a few particular values of 7 the 
function Wq{p) was calculated numerically in [|T^, around mp ~ 1 . In for mp ^ 1, 
this function was approximated by a simple expression 

This expression looks very convenient for fast estimates of various quantities. That is why 
Eq. (|T^) is widely cited and used in many papers (see, e.g., Refs. |]TT], |T^, |T3|, 0, 0, |T^, 
|2l[|). Now we show that Eq. (|T5|) is incorrect. We find out that there are two scales in 
dependence of Wo{p) on p: in the region of relatively small impact parameters, 1 -C mp < 7, 
we obtain 

Wo{p) = 1^ [2 In 72 - 3 In (mp)] In (mp) , (16) 

while in the region of relatively large impact parameters, 7 < mp <^ 7^, we have 

Note that the function Wo{p) given by Eqs. (|16|) and (|l7|) is the continuous function at 
mp = 7 together with its first derivative. Certainly, the integration of Wo{p) from Eqs. (p!6D 
and ([r?]) over p gives the main term (oc L^) in Eq.(p. 

The most important distinction between Eq. ([TB|) and our result is a quite different depen- 
dence of Wq{p) on 7 for 1 <^ mp <^ 7. It follows from Eq.(|TB|) that Wo{p) oc L at mp ~ 1 
while Eq.(|T3]) gives Wo{p) oc at mp ~ 1. Since cXunit and cr„ (for n > 2) are determined 
by the region of integration mp ~ 1, the prediction obtained with the help of Eq. ([TsD leads 
to incorrect dependence of these quantities on 7 (see below). 

Let us consider the derivation of Eqs. ([16|) and (|l3). The main contribution to the 
integrals in (|^) at p ^ 1/m comes from the region of integration 

\pj_ + q±\ ^ m , \Pz\, \qz\ ^ m-f , 

Passing to the variables P = p± + q±, r = {p^ — q^)/2, E = Pz + Qz, and x = {Pz — Iz)/ iPz + 
Qz), we take the integrals over r and x (see 0). Then, we obtain: 

9n m'^ Jm h J 
jij ^ f d^kk' {P - ky expjikp) 

J (27r)2 [fc2 + [(p _ + ^4/(^2^2)] • ^ ) 
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Wo{p) ' 



After the integration over k we come to 

m fj p fm fj p 

97r^ Jm £/ [Jl/p F Jmiix{l/p,E/y} F J 

where ^{x) is the step function. Straightforward integration leads to Eqs. (|16|), (|17|). From 
our derivation it is clear that two large logarithms at mp ^ 1 come from the integration 
over E and P, while at mp ~ 1 the logarithm from the integration over P is absent and the 
only logarithm L arises from the integration over E. 

The result (|l^), ( |T7|) can also be obtained within the standard Weizsacker- Williams ap- 
proximation. In this approximation the colliding nuclei with impact parameters pi and p2 
emit drii and dn2 equivalent photons with energies uji and u}2'- 

dui = ^ « = 1, 2 , uji,<^mj, — < Pi < — . (20) 

7r^ uji Pi m uJi 

Then these photons collide and produce a pair with the invariant mass squared s = 4ujiuj2. 
Therefore, the function Wq{p) is 

Wo{p) = J dnidn25{pi - p2 - p) ^77(5) (21) 

where o"^^ is the cross section of the process 77 e^e~ . Then we perform the integration 
over U2 using the relation 

f°° ds , , 14:71 

The main contribution to this integral is given by the region near the lower integration 
limit s = 4m?. That it is why we can extend the upper limit up to infinity and substitute 
UJ2 = m'^/uJi in the step function. After that we have 

^o{p) = Tr^— / 2f PiM — -P-Pi • 23 

Qtt'^ Jmh uJi J pi[p — piY \uJi J J 

The main contribution to this integral is given by two regions: 1/m <^ pi <^ p and 1/m ^ 

\p — Pi\ <^ p. Then we have 

28 C dui fP dpi 



OtT^ {mpy Jm/y UJi Jl/m pi 



X 



^9(7/cui -pi)^9(7CJi/m2 -p) +^9(7/cJi -p)i9(7cui/m2 -pi) . (24) 



The further integration leads again to ([16|) , ([T7|) . 

As we argued above, the function Wq{p) at mp <^ 7 and L ^ 1 has the form 

Wo{p) = CLF{mp) (25) 




FIG. 1: The function F{x) from Eq. ( pq ) vs. x = mp (solid curve) and its asymptotic form 
Fasyrap{x), Eq. (|6|)(dashed curve). 



with the universal function F{x) independent of Za.b and 7. 

We extract this function (see Fig. |1]) from the numerical results in |]I5| for the case 
7 = 3400. It follows from ( [T^ ) that the asymptotic form of the function F{x) at x ^ 1 reads 

56 In a; , . 

Using this function we obtain in the region (i) the unitarity correction to the one-pair 
production cross section 



unit 



and the cross sections for the n-pair production [n > 2) 

_ ^ 

Cr, = — F"(x)xrfx, 

n! JO 

C2 = 1-33, C3 = 0.264, C4 = 0.066, C5 = 0.0176 . 



(27) 

(28) 

(29) 
(30) 



Our result oc L" for n > 2 differs considerable from the incorrect results (T„ oc L'^^ of |TT 
and an oc L^" of 



Let us pass to the consideration of the case (ii ). If C ~ 1, but C {ZA,Baf < 1, then 
we can neglect the Coulomb corrections in W{p) but should keep the exponent in Eq.(|l^). 
It gives the result similar to Eq. (pSD for (T„ with the replacement 



27r 

Cn ^ C„(7, Za,b) = :^l F'^i^) exp[-C L F{x)] xdx . 



(31) 



For unitarity correction we have 

c 



unit 



27v-^L F(x){l-exp\-CLF(x)]} xdx. 
Jo 



(32) 



If ( {Za^bc^YL ~ 1, then we should use in Eq. (|r^) the function W{p) calculated exactly 
with respect to the parameters Z^^^a. Note that W{p), for L ^ 1 and mp ~ 1, has the 
form similar to Eq. (|25|) 

W{p) =CL F{mp, ZAa, Zbo) . (33) 

The function F{mp, Z^a, Zbo) —>■ F{'mp) for Za,_b« ~^ 0, but the difference F — F can be 
neglected in the exponential factor exp(— ly) only if this difference is small as compared to 
1 / (C-L) rather than to F. 




FIG. 2: The exact function W{x) (solid curve), Born approximation Wq{x) (dashed curve), and 
the asymptotic form of Wq{x) from (0), (0) (dotted curve) in units of C,. 



The function W{p) was calculated numerically in |^ for the particular case 7 = 100, 
Z = 79. In Fig. it is shown together with Wq{x)^ calculated also in and our 
asymptotics (0), (|T^. There is a good agreement of our analytical results with the exact 
numerical one already at mp > 2. It is also seen a noticeable difference between W{p) and 
Wq{p) in the region mp ~ 1, that is essential in calculations of the unitarity correction and 
cr„, for n > 2. Emphasize that the Coulomb corrections oc in (|lTD arise due to the 
difference between W{p) and Wq{p) in the region mp ^ 1. 



Using the numerical results for W{p) and Wq{p) from ||22[ we find that the exact value 
of <^unit /<7Born IS GQual to —4.1% , whilc the result without Coulomb effects gives —6.4% . 
With the help of Eqs.(^) and (|27|) we obtain —9.3% and —12%, respectively. Thus, in the 
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example considered (7 = 100, Z = 79) the account of the Coulomb effects as well as the 
exponential factor is very important. 

It is interesting to estimate the unitarity correction for the LHC case (7 = 3000, Z = 82). 
Using the recent numerical results of K. Hencken (private communication) for W{p) at 
7 = 3400, Z = 82, we find that the exact value of crunit/o"Born for this case is equal to —3.2%. 
This ration remains approximately the same for the LHC case. 

We are very grateful to K. Hencken for sending us the numerical data cited above. V.G.S. 
would like to thank A. Baltz, F. Gelis, L. McLerran, and A. Peshier for useful discussions 
during his stay in BNL. This work was supported through Grants RFBR 00-02-17592, 01- 
02-16926, and St. -Petersburg E 00-3.3-146. 
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